NRG Study of an Inversion-Symmetric Interacting Model: 
Universal Aspects of its Quantum Conductance 
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We consider scattering of spinless fermions by an inversion-symmetric interacting model char- 
acterized by three parameters (interaction U, internal hopping td and coupling t c ). Mapping this 
spinless model onto an Anderson model with Zeeman field, we use the numerical renormalization 
group for studying the particle-hole symmetric case. We show that the zero temperature limit is 
characterized by a line of free-fermion fixed points and a scale r(U,t c ) of td for which there is perfect 
transmission. The quantum conductance and the low energy excitations of the model are given by 
universal functions of td/r if td < T and of td/t^ if td > I\ T = tj. being the level width of the 
scatterer. This universal regime becomes non-perturbative when U exceeds F. 

PACS numbers: 71.10.-w,72.10.-d,73.23.-b 
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In quantum transport theory, the conductance G of 
a nanosystem inside which the electrons do not interact 
is given by g — G/(e 2 /h) = \t ns \ 2 when the tempera- 
ture T — > 0, \t ns \ 2 being the probability for an electron 
at the Fermi energy Ef to be transmitted through the 
nanosystem. This Landauer-Buttiker formula can be ex- 
tended to an interacting nanosystem, if it behaves as a 
non-interacting nanosystem with renormalized parame- 
ters. We study such a renormalization using the numer- 
ical renormalization group (NRG) algorithm [J, Q and 
an inversion-symmetric interacting model (ISIM) which 
describes the scattering of spin-polarized electrons (spin- 
less fermions) by an interacting region characterized by 
an internal hopping term td, a coupling term i„ and an 
interaction strength U. This model was used [3, 0] for 
studying the effect of an external scatterer upon the effec- 
tive transmission of an interacting region, assuming the 
Hartree-Fock (HF) approximation. We revisit ISIM with 
the NRG algorithm for investigating non-perturbative 
regimes where other methods (NRG or DMRG algo- 
rithms) than the HF approach become necessary. 

Quantum impurity models [H, as the Anderson model 
which describes a level with Hubbard interaction U cou- 
pled to a 3d bath of free electrons, were introduced to 
study the resistance minimum observed in metals with 
magnetic impurities. The Kondo problem refers to the 
failure of perturbative techniques to describe this min- 
imum. The solution of these models by the NRG al- 
gorithm, a non-perturbative technique 0, 0] introduced 
by Wilson, is at the origin of the discovery of univer- 
sal behaviors which can emerge from many-body effects. 
The observation [f| of the Kondo effect in semiconduc- 
tor quantum dots has opened a second era for quantum 
impurity models, now used for modeling mesoscopic ob- 
jects (single [1} or double Q quantum dot systems) in- 
side which electrons interact, in contact with baths of 
free electrons (large conducting non interacting leads). 

Though the Kondo effect is induced by magnetic mo- 
ments, it is also at the origin of spinless models, such 



as the interacting resonant level model [8[ (IRLM) which 
describes a resonant level (Vdcftd) coupled to two baths of 
spinless electrons via tunneling junctions and an interac- 
tion U between the level and the baths. IRLM, which is 
often used for studying nonequilibrium transport [1, @] , is 
related to the Kondo model, the charge states rid = 0, 1 
playing the role of spin states. Both ISIM and IRLM 
are inversion symmetric and can exhibit orbital Kondo 
effects. However, the Zeeman field acting on the impu- 
rity is played by the hopping term td for ISIM, and by 
the site energy Vd for IRLM. Therefore, ISIM does not 
transmit the electrons without field, while IRLM does. 
The two-particle states have been given for ISIM (loj . 

For the particle-hole symmetric case Q , the Anderson 
model maps onto the Kondo Hamiltonian if U > irT, T 
being the impurity-level width. In that case, there is a 
non-perturbative regime where the temperature depen- 
dence of physical observables such as the impurity sus- 
ceptibility is given by universal functions of T/Tk, Tk 
being the Kondo temperature. If U < 7iT, the impurity 
susceptibility can be obtained by perturbation theory. 
Mapping ISIM onto an Anderson model with a Zeeman 
field td, and assuming that the role of td should quali- 
tatively resemble that of a finite temperature, we expect 
the following scenario for the ISIM conductance g of the 
particle- hole symmetric case: If U > ttT oc t 2 , we expect 
a non-perturbative regime where g should be given by 
a universal function of td/r independently of the values 
of U and t c , with a scale r(t c , U) of td playing the role 
of a Kondo temperature Tk- If U < ttT, the HF theory 
should correctly give g. This scenario will be more or less 
confirmed by extensive NRG calculations. 

ISIM Hamiltonian: H = H ns + Hi + H c . The Hamil- 
tonian of the interacting region (the nanosystem) reads: 



td (cjq + c\c^j + Vg {no + m) + Un ni . (1) 



Hns — 



cj. and c x are spinless fermion operators at site x and 
n x = c x c x . The leads are described by an Hamiltonian 
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FTC!. 1: Line of free fermion fixed points ((7 = 0, thick solid 
line) characterizing ISIM when T — > as td increases from 
td = (SC fixed point) towards td — > oo (PO fixed point). 
The FO, LM and SC fixed points and the RG trajectories [l[ 
followed by ISIM as T decreases for td = are indicated in 
the plane td = 0, for 7rr > U (dashed) and 7rr < U (solid). 

E/OO / 
( c x c x+i + H.c), where } means 

that x = — 1, 0, 1 are omitted from the summation. The 
coupling Hamiltonian H c = — t c (c_iCQ + c\c 2 + H.c ). 

Mapping onto an Anderson model with Zeeman field: 
Because of inversion symmetry, one can map ISIM onto 
a semi-infinite Id lattice where the fermions have a 
pseudo-spin and the double site nanosystem becomes a 
single site with Hubbard repulsion U at the end point 
of the semi-infinite lattice. %i ox = ( c _x+i ± c J)/v^ 
creating a spinless fermion in an even/odd (e/o) com- 
bination of the orbitals at the sites x and — x + 1 
of the infinite lattice, (or a fermion with pseudo-spin 
(j = e/o in a semi-infinite lattice), one gets H ns — 
(Vg - t d )n c + (Vg + t d )n + Un c n a , where n a = al 1 a a l 
and where the pseudo-spin "e" ("o") is parallel (anti- 
parallel) to the "Zeeman field" td- In terms of the 
operators d\ = J2x°=2 sin(fc(x — l))aj. x creat- 

ing a spinless fermion of pseudo-spin a and momen- 
tum k in the semi-infinite lattice, Hi = ^2 k £knt,a and 

H c = J^k a V(k)( a t l^fe o-+-^- c -)i where the fc-dependent 
hybridization V(k) = —t c y/2/w sink yields an impurity 
level width T = t 2 , nk,a — d\ a d k and e/. = — 2t^cosfc. 
ISIM is almost the Anderson model, except that the im- 
purity has a Zeeman field td and is coupled to a semi- 
infinite Id bath of free electrons. When t d — > 0, ISIM 
exhibits an orbital Kondo effect if the equivalent Ander- 
son model can be reduced to a Kondo model. 

NRG procedure: ISIM can be studied using Wilson's 
procedure 0, Q developed for the Anderson model after 
minor changes. First, we assume V(k) as V(kp = tt/2) 
and, taking A = 2, we divide the conduction band (log- 
arithmic discretization) of the electron bath into sub- 
bands characterized by an index n and an energy width 
d n = A~ n (l — A -1 ). Within each sub-band, we introduce 
a complete set of orthonormal functions ipnp{t), and ex- 
pand the lead operators in this basis. Dropping the terms 



with p / and using a Gram-Schmidt procedure, the 
original Id leads give rise to another semi-infinite chain 
with nearest neighbor hopping terms, each site being la- 
belled by the same index n as the energy sub-band from 
which it comes, and representing a conduction electron 
excitation at a length scale A™/ 2 ^ 1 centered on the im- 
purity. In this transformed Id model, the successive sites 
are coupled by hopping terms t n ^ n+ \ cx A-™/ 2 which van- 
ish as n — > oo. The impurity and the N—l first sites form 
a NRG chain of length N and of Hamiltonian H^. This 
length can be interpreted Q as a logarithmic temper- 
ature scale. The NRG chain coupled to the impurity is 
iteratively diagonalized and rescaled, the spectrum being 
truncated to the N s first states at each iteration. The be- 
havior of ISIM as T decreases can be obtained from the 
spectrum of Hn as N increases, the bandwidth of Hn 
being suitably rescaled at each step. A fixed point of 
the RG flow corresponds to an interval of successive even 
(or odd) values of N where the rescaled many-body ex- 
citations Ei(N) do not vary. If it is a free- fermion fixed 
point, Ei = J2 a e a, the e a being one-body excitations, 
and the interacting system behaves as a non-interacting 
system (U — 0) with renormalized parameters td and t c 
near the fixed point. Moreover, if one has free fermions 
when T — > 0, g can be extracted from the NRG spectrum. 

Symmetric case: Using this NRG procedure, ISIM can 
be studied as a function of T for arbitrary values of its 
bare parameters. Hereafter, we take th = 1, Ep = 
and Vg = —U/2. This choice makes ISIM invariant 
under particle-hole symmetry, with a uniform density 
{( n x) = 1/2) and 3 effective parameters (U,t c ,t d ). 

Suppression of the LM fixed point as td increases: 
When td — 0, ISIM is an Anderson model which has 
the RG flow sketched in Fig. [1] for the particle-hole sym- 
metric case. At low values of N (high values of T), ISIM 
is located in the vicinity of the unstable free orbital (FO) 
fixed point. As N increases (T decreases), ISIM flows 
towards the stable strong coupling (SC) fixed point. If 
7rf 2 < U, the flow can visit an intermediate unstable fixed 
point: the local moment (LM) fixed point before reach- 
ing the SC fixed point. In that case, ISIM is identical to 
a Kondo model characterized by a temperature Tk and 
by universal functions of the ratio T/Tk- If 7ri 2 > U, 
the flow goes directly from the FO fixed point towards 
the SC fixed point, and there is no orbital Kondo effect 
for t d — > 0. In Fig. [D^a), the first many-body excitations 
Ei of ISIM are given for increasing even values of N for 
td = 0. Since 7rf 2 < U, one gets 3 plateaus correspond- 
ing to the 3 expected fixed points. Inside the plateaus, 
the spectra are free-fcrmions spectra which are described 
in Ref. 0]. However, between the plateaus, there are no 
free- fermion spectra and Ei ^ ^2 a e a- As td increases 
(Fig. [2 a)), the LM plateau decreases and vanishes when 
td~U, 

Evolution of the SC fixed point as td increases: In 
the limit N — > oo (T — > 0), let us study the Ei as a 
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FIG. 2: (Color online) Fig. [2ja): Many body excitations _B/ as a function of N (even values) for U = 0.005 and t c = 0.01. 
For td — (Fig. [2Jal)), one can see the 3 successive plateaus (FO, LM and SC fixed points) of the Anderson model 2|. As td 
increases (Fig. [2ja2) and Fig. EJa3)), the LM plateau shrinks and disappears when td ~ U. Fig. EJb): One body excitations 
£ a (td) (extracted from the Ei(N — > oo, td)) for [/ = 0.1 and t c = 0.1 (left scale). The solid (dashed) line corresponds to NRG 
chains of even (odd) length N. Conductance g(td) extracted from Ae(td) using Eq. ([2]) (thick red curve, right scale). For td = T, 
the e a are independent of the parity of N and g = 1. Fig. [3Jc): For [7 = 0, e a (td/tc) an d g(td/t^) extracted from the NRG 
spectra (x). g — cosh _2 (X) (red line) with X = ln(td/tc) is correctly reproduced. Figs. [2jd),(e): g{td) for t c = 0.1 and many 
values of U, calculated by NRG algorithm (d) and by HF theory (e). In Fig. EUd), the larger is U, the smaller is td = r where 
3 = 1. The curves correspond respectively to U = 0.25, 0.2, 0.15 (3 left peaks) and U = 0.1, 0.09, . . . , 0.01, (11 right peaks). 
In Fig. [2je), the HF values are accurate for U = 0.02, 0.01, (3 right peaks), but become inaccurate when U w 0.04 w T. For 
U > F, the HF curves (dashed lines) are very different of the corresponding NRG curves (Fig. EJd)). 



function of t d . For td = 0, one has the SC limit [2] 
where the impurity is strongly coupled to the second site 
(the conduction-electron state at the impurity site) of 
the NRG chain. The impurity and this site form a sys- 
tem which can be reduced to its ground state (a sin- 
glet), the N — 2 other sites carrying free fermions exci- 
tations e a which are independent of that system. In the 
presence of a Zeeman field td ^ 0, the free-fermion rule 
Ei{td) = J2 a £a{td) remains valid (see Fig. [DJb)) and 
the T — > limit of ISIM is given by a continuum line 
of free-fermion fixed points where U — 0, as sketched in 
Fig. [TJ When the pseudo-spin degeneracy is broken, the 
first (second) one-body excitation ei (£2) carry respec- 
tively an even (odd) pseudo-spin if TV is even. This is the 
inverse if N is odd, e\ (52) carrying respectively an odd 
(even) pseudo-spin. For td — > 00, the impurity occupa- 
tion numbers n G = 1 and n Q = 0, and the N — 1 other 
sites of the NRG chain are independent of the impurity. 
We call this fixed point "Polarized Orbital" (PO), since it 
coincides with the FO fixed point of the Anderson model, 
except that the spin of the free orbital is fully polarized 
in our case. Since for N — > 00 and t d — > (SC fixed 
point), the free part of the NRG chain has N — 2 sites, 
while it has N — 1 sites for td — > 00 (PO fixed point), 
there is a permutation of the e a (td) as td increases: as 



shown in Figs. H{b) and (c), the e Q (trf — ► 0) for even 
become the e a {td °°) f° r N odd and vice- versa. 

Characteristic energy scale t: We define the charac- 
teristic energy scale r(t c ,U) of ISIM as the value of td 
for which the e a (td) are independent of the parity of N 
when N —> 00. Because of particle-hole symmetry, the 
nanosystem (the impurity of the NRG chain) is always 
occupied by one electron. Binding one electron of the 
leads with this electron reduces the energy when td < T, 
while it increases the energy when td > t. For td = t, it 
is indifferent to bind or not an electron of the lead with 
the one of the nanosystem, making ISIM perfectly trans- 
parent. This gives the proof that, for every values of U 
and t c , there is always a value r of td for which g = 1. 
The argument is reminiscent to that giving the condition 
for having a perfectly transparent quantum dot in the 
Coulomb blockade regime: td in our case, the gate volt- 
age in the other case, have to be adjusted to values for 
which it costs the same energy to put an extra electron 
outside or inside the dot. 

Extraction of the conductance g from the NRG spectra: 
If <5 e (S ) are the even (odd) scattering phase shifts at Ep, 



g(t d ) = sin 2 (<5 e - 5 a ) 
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FIG. 3: (Color online) Conductance g as a function of td/r for 
3 values of t c and many values < U < 35. Inset: T(U,t c )/U 
as a function of U /t 2 (+) and fit r = t 2 exp —(U /(7rf 2 )) (solid 
red line). 

where Ae = ti — 1\ is the energy gap between the two first 
excitations of a NRG chain of even length N — > oo (see 
Fig. E](b) ) . When U — 0, this relation is a consequence of 
Friedel sum rule, which can be written for each pseudo- 
spin channel separately. In that case, g = COSll (A) 
where X = \n{td/t 2 c ) and r = t%. The e a (td) given by the 
NRG algorithm for U = are shown in Fig. [2jc) with the 
corresponding values of g obtained from Eq. ((2|) , showing 
that this proceduregives correctly g when U = 0. It has 
been shown BlIIEl thatEq. © can also be used when 
U 7^ 0, if there are free fermions when T — > 0. 

Non-perturbative regime (U > T/A): In HF theory, td 
takes [3] a value v = td + J7(cqC 1 (u, t c )) and g = 1 if 
v = t 2 c . This gives for the scale r a HF value thf — 
t 2 — AU where A = (cjc^t; = t 2 ,t c )) depends weakly 
on t c , A = 1/tt (1/4) for t c = 1 (0). When U -> t 2 JA, 
thf — ► 0, showing that HF theory cannot be used above 
an interaction threshold which is almost the threshold 
7rF giving the onset of the non-perturbative regime for 
the Anderson model. This breakdown of HF theory for 
U w T/A can be seen if one compares Fig. [2jd) (NRG 
results) and Fig. \2le) (HF results). 

Universality: The conductance g extracted from the 
NRG spectra for t c = 0.01,0.1 and 1 and < U < 35 
is given as a function of td/r in Fig. [3] One can see 3 
successive regimes. When td < r, there is a single curve 
which is independent of U and t c and which corresponds 
to g = cosh~ 2 (JT) with X = ln(td/V), and not \n(td/t 2 ) 
as for U — 0. When td > t, another universal curve 
independent of t c and U describes the data as a function 
of td/r as far as td does not exceed T. Indeed, the same 
data plotted as a function of td show that g becomes 
independent of U when td > T. In this third regime 
(parallel lines which can be seen in Fig. [3] for large values 
of td/r) g = coslT 2 pO with X = \n(t d /t 2 c ) as if U = 0. 

Roles ofT and td: We have assumed analogies between 
the effect of T in the Anderson model, the effect of a Zee- 



man field at T = in the Anderson model, and eventually 
the effect of td at T = in ISIM. This was based on the 
idea that the singlet state of the SC limit could be broken 
either if the temperature T or the Zeeman energy td ex- 
ceeds T/f. Let us summarize the interest and the limit of 
these analogies. Increasing T in the Anderson model (or 
in ISIM with td = 0), one gets 3 regimes, each of them 
being characterized by a single fixed point (Fig. [2{a)). 
There are no free fermions for temperatures T » Tk 
(SC-LM crossover) and T s» T (LM-FO crossover). In 
contrast, increasing td in ISIM at T = 0, one has al- 
ways free fermions (Fig. [2b)), and not only around 3 
fixed points. However, there are 3 regimes in ISIM as 
td increases, as in the Anderson model as T increases, 
delimited by 2 energy scales r and T. The behavior of 
t w t 2 exp — (U/{-Kt 2 )) (inset of Fig. [3]) resembles that 
oiT K fa i cV ^r!772exp-(C//(7rt 2 )) (in ISIM units), while 
the second scale is given by T in the 2 models. Eventu- 
ally, we point out the similarity between the universality 
discussed in this letter for g and that which character- 
izes [3| also at T = the behavior of the singlet-triplet 
gap for a magnetic impurity confined in a box of mean 
level spacing A, as a function Tk/A- 
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